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We analyze a simple Split Supersymmetry scenario where fermion masses come from 
anomaly mediation, yielding m s ~ 1000 TeV, ~ 100 TeV, and nrif ~ 1 TeV. We 
consider non-thermal dark matter production in the presence of moduli, and we find 
that the decay chains of moduli — > LSPs and moduli —* gravitinos —* LSPs generate 
dark matter more efficiently than perturbative freeze-out, allowing for a light, LHC 
visible spectrum. These decaying moduli can also weaken cosmological constraints 
| on the axion decay constant. With squark masses of order 1000 TeV, LHC gluinos 

■ will decay millimeters from their primary vertices, resulting in a striking experimental 

signature, and the suppression of Flavor Changing Neutral Currents is almost sufficient 
to allow arbitrary mixing in squark mass matrices. 
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p_; 1 Introduction 

<D '. 

With the Large Hadron Collider only a few years away, a new division has arisen in the particle 
physics community - is the unnaturalness of the standard model a problem to be solved 
through clever model building, or is it a hint that physics at the TeV scale is different than 
we have imagined, so that the values of dimensionful parameters are determined by anthropic 
fine tuning? As physicists, we would prefer a 'physical' explanation of small numbers such 
as the cosmological constant and the higgs mass, but the success of Weinberg's prediction of 
the cosmological constant p] and the level of fine tuning necessary even in our best theories 
forces us to take the anthropic argument seriously. We will be taking it seriously for the 
present work. 

Yet the mind of an anthropically-motivated model builder is a troubled one, for it seems 
that this profound shift in our worldview has had the unfortunate side effect of putting our 
whole enterprise out of business. If the weak scale is determined by anthropic selection, 
perhaps there is no new physics at the TeV scale, and the LHC and even the ILC will be 
colossal disappointments. 

There are two simple reasons to remain hopeful about new TeV scale physics: gauge 
coupling unification and dark matter. The former requires new particles charged under the 
standard model gauge group, and the latter requires a new stable particle with the correct 
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relic density. If we assume for maximal simplicity that these two problems have a common 
solution, then it is reasonable to assume that our dark matter candidate will have a weak 
interaction cross section. Furthermore, if we assume that the relic dark matter density is 
determined by perturbative freeze-out, then we can expect new physics at the TeV scale, 
tightly constrained by gauge coupling unification. Some of these assumptions may be wrong, 
but that is not necessarily a bad thing - the minimal model that follows from this reasoning 
will probably be invisible to the LHC [22J. 

A very popular set of models that result from this methodology fall under the heading 
of Split Supersymmetry [2], where the scalar superpartners necessitated by supersymmetry 
are very heavy, while the fermions, protected by chiral symmetry, lie near the TeV scale. 
The purpose of the present work is to show that in one particularly simple model of Split 
Supersymmetry, perturbative freeze-out is not the dominant mechanism for generating dark 
matter. However, the mechanism that will replace it is more efficient, allowing for an even 
lighter, more LHC-visible spectrum. This seems like good news, because in standard Split 
Supersymmetry scenarios [3], [IS], dark matter is expected to be too heavy to be seen at 
the LHC; however we will see that in our case there is a danger of overclosure. One of the 
most elegant scenarios is ruled out, but we explore several ideas that can rescue it, and the 
scenario of Moroi and Randall [13] remains a viable possibility. 

To understand our mechanism, we first need to explain the spectrum of the model. As 
shown in [6], [19], wherever supersymmetry is broken, there will be visible sector supersym- 
metry breaking effects from Anomaly Mediation. Furthermore, the methodology of effective 
field theory requires that we include in the Lagrangian all higher dimension operators al- 
lowed by symmetry, with appropriate suppression by inverse powers of the cutoff. When 
supersymmetry is broken in a hidden sector, F and D term VEVs in these operators will in- 
duce supersymmetry breaking in the visible sector. These two effects are completely generic, 
requiring no theoretical gymnastics, so a model where supersymmetry breaking arises only 
in this way would be particularly elegant. This is the model that we will analyze; it was 
studied for related reasons in [18J. 

With Anomaly Mediated Supersymmetry Breaking, gaugino masses are given by a loop 
suppression factor times m 3 / 2 , so we expect that m 3 / 2 ~ 100 TeV so that new fermions are 
near a TeV. Contributions from higher dimension operators suggest that the scalar mass scale 
m s ~ 1000 TeV. In supersymmetric theories, there generically exist weakly interacting moduli 
that get masses of order m s after supersymmetry breaking. Now we see the complexity of 
the cosmology - not only do we have to account for the perturbative freeze out of the LSPs, 
but we also need to consider processes involving late-decaying moduli and gravitinos. These 
late decaying particles produce a great deal of entropy, potentially weakening cosmological 
constraints on the properties of axions. 

We will see that with our mass spectrum, dark matter production is dominated by two 
potential decay chains: moduli — ► LSPs and moduli — > gravitinos — > LSPs. Only the first 
process is available for < 2m 3/ / 2 J^| whereas both decay chains are open for > 2m 3 / 2 - 
thus the physics is qualitatively different for different values of the modulus mass. The first 
process was studied in [13], and it naturally leads to the correct relic density in our scenario 
when rrifp < m 3 / 2 . The most naive version of the second decay chain is too efficient [33], [31], 

2 Of course this depends on the available gravitino- producing modulus decay channels - the strict limit is 
777,0 < 7773/2 ~~ but under reasonable assumptions <f> — > 2ip 3 /2 is the only available decay mode. 
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[35] , so we consider a variety of mechanisms to alleviate this problem, including symmetries in 
the modulus sector, scenarios with multiple moduli, KKLT type moduli, and supersymmetry 
breaking dominated by D-terms. Of course we do not need these mechanisms for a successful 
cosmology if < ^3/2, but it is interesting to explore all reasonable options. 

Our model also has an exciting LHC signature. The gluino must decay through a virtual 
squark, so with 1000 TeV squarks, it is very likely that the LHC will see displaced gluino 
vertices if gluino production is kinematically allowed. In fact, we expect that gluinos will be 
copiously produced, since cosmology suggests a light spectrum. 

We are also in an interesting region for Flavor Changing Neutral Currents. The strongest 
constraints from FCNC come from the e parameter of the K-K system, which depends on 
the imaginary parts of the mass insertion parameters. Setting all mass insertions equal to a 
common value S (there are no significant cancellations), we find from [7] that 

Sr /1000 TeV\ 2 
\ m s J 

Experiments constrain e < 2 x 10~ 3 , so for m s pa 1000 TeV, we find that as long as 5 < 1/2, 
our scalar mass scale does not conflict with observations of FCNC. We also expect that electric 
dipole moments induced by new interactions may be visible at next generation experiments, 
as shown in [3], and our spectrum may help to explain neutrino masses, as examined in [18] . 

The outline of this paper is as follows. In section two we display the mass spectrum of 
our model. In section three we consider the cosmological implications of moduli, gravitinos, 
and LSP (Wino) dark matter. In section four we show that displaced gluino vertices will be 
a generic LHC signature for our model, and with section five we conclude. 



2 The Mass Spectrum 

We obtain the fermion mass spectrum 

m 3/2 = M = (F > ~ 100 TeV (2) 

m- g = ^(F,)~lTeV (3) 

from anomaly mediated supersymmetry breaking [6], [19] . The details of hidden sector su- 
persymmetry breaking are unimportant because the hidden sector only communicates with 
standard model particles through the auxilary field of supergravity and through higher di- 
mension operators in the Kahler potential. We assume that the scalars acquire masses from 
operators such as 

L D / d^—^V WQ (4) 

J M GUT 

where X is a hidden sector field, because these terms cannot be forbidden by any symmetry. 
We expect these operators to arise generically when we integrate out GUT or string scale 
particles, making them a factor of 10 - 100 times larger than 7723/2. Thus we can estimate 
that 

m s = m 3/2 ^- ~ 1000 TeV, (5) 



GUT 
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completing our rough picture of the mass spectrum. 

It is also important to consider the generation and the effects of the /i and Bfi terms. 
Will n and Bfi be near m s , WI3/2, or at the TeV scale? First, there can be contributions from 
R-symmetry breaking spurionic operators X = 1 + 9 2 m s such as 

L d J d i 9eX^XH 1 H 2 (6) 

L D J d 4 9eX j H 1 H 2 (7) 

where the factor of e is included to parameterize a possible approximate PQ symmetry. 
This gives rise to B[i ~ em 2 and fi ~ em s . In most split supersymmetry scenarios we 
use R-symmetry to prohibit such spurionic contributions because they produce a nearly 
degenerate spectrum. However, in our scenario m^n/mf ~ 100, which can be conveniently 
explained as a loop factor from anomaly mediation. Thus such operators may be permitted. 
If supersymmetry breaking does not lead to R-symmetry breaking we only have spurionr^ 
such as 

L D / d 4 9eYHiH 2 (i 



with Y — 1 + 9 4 m 2 , which contribute to Bfi only. Finally, there are contributions from the 
conformal compensator of supergravity 

L D J d 4 0e0V#i#2 (9) 

where = l + 9 2 m^/2- l so that [i ~ Thus we find that Bfi pa em 2 in all cases, and either 

\l pa em s or fi pa . We can take e ~ 1/100 to explain the top-bottom mass hierarchy, 

giving a /i term at the TeV scale. 

The presence of the \x term modifies the gaugino masses [6], [19j so that 
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m 3/2 (10) 

13/2 (11) 

(12) 



where m^sin(2/?) = Bfi and the function 



/M = ^ (13) 

X — 1 

For small values of 2 (the expected case) this is a small effect. It increases the separation 
between the wino, bino, and gluino mass scales, but there is no real qualitative change in 
the spectrum. For large x we can obtain a very light Wino LSP, although we would need to 
abandon the PQ symmetry; this will be breifly considered in section 3.3.3. 



3 As explained in [3], this does not require that supersymmetry breaking comes from gauge superfields, 
but only that supersymmetry breaking is not accompanied by R-symmetry breaking. 
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3 Cosmology 



Our variant of the Split Supersymmetry mass spectrum includes 100-10, 000 TeV moduli 
field(s) and a 100 TeV gravitino, so we must check that these new ingredients do not 
disturb Big Bang Nucleosynthesis, and that the correct relic abundance of LSP dark matter 
candidates obtains. We begin with some general issues about our universe's history, and then 
we consider the possibility that our moduli fields weaken cosmological bounds on the axion. 
Finally, we explain the details of our proposed LSP generation mechanism. 

We can think of our modulus field as though it were an inflaton - it begins with a VEV 
0o that can be of order M p i. Once the hubble constant decreases to H ~ m^, the modulus 
begins to oscillate and its energy density red-shifts like that of matter. This occurs very 
early in the universe's history, with T ro u ~ w 1 M p im^ ~ 10 9 TeV. At a time t eg shortly after 
4> begins to roll (assuming <f> begins with a Planckian VEV), p^ = pR a d, and from this point 
until t decay = Tt ] the universe is modulus dominated. We take the modulus decay rate as^| 

r« = —4 = (10- 5 sec)" 1 (N) ( m * Y (14) 
16ttM p 2 z V ' \10J \ 1000 TeV/ V ; 

where M v \ = 2.4 x 10 18 GeV is the reduced Planck mass. We can think of N as the number of 
available light decay modes, although it is strongly dependent on unknown, 0(1) couplings. 
We see that the modulus decays well before BBN. The decaying modulus will dilute particles 
left over from the original inflaton decay, reheat the universe, and re-populate the universe 
with gravitinos. 

Any thermal gravitinos produced during a prior inflationary reheating are diluted by the 
entropy from the modulus decay, which is given by 



3/4 



M, 



$after ( P<f>(t decay) \ 

Sbefore \PR(tdecay) J TTlji \M] 



(15) 



where the ratio of the energy densities at the time of modulus decay simply comes from the 
two different equations of state. If we assume 0o is m the neighborhood of M p i, then the 
entropy released is very large. Modulus decay reheats the universe to a temperature [UJ 

T^l.l^yr^=180MeVx ( ^ f^-Y^f ^ f 2 (16) 
H y* V tot P i y 1Q j Vl0.75/ VIOOOTeV/ V ; 

which is very low, so we certainly do not thermally regenerate heavy particles. To avoid 
disturbing BBN, we need Tr ^> 1 MeV, so 3> 30 TeV. Similarly, the gravitino lifetime is 

9 /100TeV\ 3 . , 

r 3/2 = 5 x 10~ 2 sec 17 

V m 3/2 J 



4 There is certainly a significant uncertainty because we do not know any of the 0(1) couplings involved. 
However, we do know that modulus decays to matter in the complex representation of a symmetry group 
[13], [33], [35] receive a suppression proportional to a power of m mat /m^, which means that moduli decay 
dominantly to gauge bosons, gauginos, and possibly the heavy scalars. 
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so for > 60 TeV gravitinos do not disturb BBN |12j . [33] even if they are re-introduced 
by modulus decay. Note, however, that the lifetime is increased in split supersymmetry 
because the gravitino cannot decay to states involving the heavy scalars. 

Dark matter is produced in two different ways: directly from decaying moduli, and indi- 
rectly from decaying gravitinos. Thermal relic gravitinos are diluted to negligible levels by 
modulus decay, so gravitinos can only arise from the decay of heavy moduli with > 2m^/2- 
Thus there are two qualitatively different types of moduli - those with > 2m 3 / 2 and those 
with rritf, < 2m^/2- In the following subsection we consider direct modulus-LSP interactions, 
focusing on moduli with < 2ms/2, and commenting on heavier moduli. Then in the next 
subsection we consider gravitino production from moduli with > 2m 3 / 2 . It turns out 
that dark matter production from gravitinos is more efficient than direct production from 
modulus decay; this is essentially because the dark matter from gravitinos is produced later 
in the universe's history, so it accounts for a larger fraction of the universe's energy density. 



3.1 Direct Modulus Decay to Dark Matter 

We expect that the modulus will have generic string or Planck suppressed couplings to MSSM 
fields. For example, we may have interactions such as 



d 2 6-^-W Q W a L d 4rFlv + ^-AA, 

M pl M pl M pl 



where A may be any gauginc@. If $ has a supersymmetric mass, then cx m^0, and both 
terms give large decay rates, but in other situations where only arises after supersymmetry 
breaking, we may have < and the fermionic decay mode will be suppressed by the 
fermion mass. Due to R-Parity conservation, this implies that an order one fraction of 
the 'modulus particles' may eventually produce an LSP. The number density of 'modulus 
particles' relative to the entropy just after reheating is roughly 

\ -1 / \ 1/2 



Y *=W ~W (19) 

If this were the whole story, we would overdose the universe. If the branching ratio to 
fermions is suppressed by mj/m^ an acceptable dark matter density might be possible. 

However, the LSPs produced by modulus decay pair annihilate [13] until their annihilation 
rate is less than H, after which they freeze out. Although the Wino's produced by the 
decaying modulus are ultra-relativistic, they acquire a thermal distribution with temperature 
Tr from interactions with the radiation bath. For instance, the process W°v — > W°u proceeds 
more quickly than H for Tr ^> 1 MeV, rendering the Winos non-relativistic. The annihilation 
cross section for Wino dark matter [13] into W ± pairs is 



, , 9 2 1 (1 - xwfft 

= 2^l2^)F (2°) 



5 Note that if $ acquires a large F-component VEV, then this operator could give dangerously large gaugino 
masses, so we must assume that either F$ is small or that the coefficient of this operator is suppressed. 
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in the non-relativistic limit, with xw = m wl m2 LSP- Thus the maximal dark matter relic 
density is reduced to 

3H 



n L SP 



2(va) 



T 2 

R (21) 



Tr 



M pl (va)' 



Note that this will apply to any late decay that produces dark matter, including gravitino 
decay. Thus we find [13] 

P - ~ 4 x 10- 10 GeV x (-] r^r 1/4 {2 ~ XW? ( mLSP Y r i00TeV V /2 (22) 

where we have taken g* = 10.75. We need p/s m 4 x 10~ 10 GeV [36] to account for dark 
matter, so if there is a single light modulus with < 27B3/2, we can easily obtain the correct 
relic abundance [13J. Since the decay of this modulus will not produce gravitinos, we have a 
viable cosmology. 

If ^> then after modulus decay we are left with a negligable density of LSPs, 
so next we consider the consequences of gravitino production in scenarios with such a heavy 
modulus. 



3.2 Moduli-Gravitino Interactions and Indirect Dark Matter Pro- 
duction 

There are two potential channels in which modulus decay produces gravitinos, 

<P -> + ^3/2 (23) 

-> 2^/ 2 (24) 

where <fi is the modulino, the superpartner of the modulus field. We do not expect the 
modulus to decay to a gravitino and another fermion because the other fermion would reside 
in the hidden sector, and we expect that hidden sector fields are very massive. For now we 
will assume that the modulino is very heavy, so that the first channel is also forbidden by 
kinematics - later we will see that relaxing this assumption does not help us to obtain a 
viable cosmology. In this section we consider the limit that ^> 7713/2, and we begin with 
the two gravitino channel. 

Gravitinos interact primarily through their longitudinal (helicity ±1/2) components, 
which correspond to the goldstino, the field that is 'eaten' in the superhiggs mechanism. 
Thus the goldstone boson equivalence theorem becomes the goldstino equivalence theorem, 
which simply tells us that amplitudes involving longitudinal gravitinos can be computed us- 
ing the goldstino instead. In the appendix it is shown that the only two derivative operator 
involving a scalar and two non-linearly realized, derivatively coupled goldstinos (longitudinal 
gravitinos) is 

l^d^x^d^x ~ jdvs^a^d^chi). (25) 
where the 'sgoldstino' field is 
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and F 2 = Eji^ 2 + £ a -D 2 gives the goldstino decay constant. The sgoldstino is generally not a 
mass eigenstate - it is a blend of many different scalar fields. Thus the rate for moduli decay 
to goldstino pairs is determined entirely by the overlap of the mass eigenstate modulus with 
the sgoldstino. 

In supergravity one expects higher dimension operators to be present in the Kahler po- 
tential, with suppression of order the GUT, string, or Planck scale. These operators tend to 
shift the VEVs of scalar fields and auxiliary fields. We would expect the effect on a modulus 
to be proportional to l/m^ because the mass tends to stabilize the field, so on dimensional 
grounds we can estimate that 

F<s> „ J^L. „ 0*11 F (27) 
m^Mpi 

for the auxiliary field of a heavy modulus. A more rigorous argument for this conclusion was 
given in [341, and it seems to be generically true in specific models in the absence of special 
symmetries^ Assuming that the cf) field really is a mass eigenstate, this implies a rate for 

13/2 = i^mJ (28) 

where C is a dimensionless constant. The branching fraction for modulus decay to gravitino 
pairs is therefore C/N. 

In the approximation of instantaneous modulus decay, the ratio of modulus particles 
before the decay to entropy immediately after the decay is simply = T^jra^. Thus we 
estimate that 



Nm^ y/N V 1000 TeV, 

Now conservation of R-parity implies that each gravitino will eventually decay into at least 
one LSP, so we can approximate Y^sp = Y$/i- If the LSPs do not pair annihilate, then we 
find 

PLSP _ R „ 1n - 6 n ^r <±_ ( m LSP \ ( ™* A 1/2 



6xlO- fa GeV^ ^ (30) 

~ V100 GeV/ V1000 TeV/ 



so we must take C/yN < 10 -4 to avoid overdosing the universe. This seems to be a very 
stringent constraint - it rules out the most elegant scenario with C ~ 1 - so in the next 
section we will consider whether a small C can be achieved naturally, or if other processes 
can reduce the relic density of LSPs. Note that the relic density is proportional to misp, so 
if this scenario can be made viable, then we expect a light spectrum. 

One might hope that the LSPs produced by gravitino decay pair annihilate, decreasing 
their abundance to within acceptable levels. We can use equation (121"]) together with the 
spectrum given in equation (1101) to write the resulting LSP density in terms of m 3 / 2 alone 



6 It continues to be true when the modulus has a shift symmetry, although as we will discuss later there 
are exceptions in some models, such as specific versions of KKLT |35j . 
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and this would require 7773/2 < 20 TeV to obtain the correct relic density. Such a light 
gravitino would cause severe problems with BBN, and a Wino with rriLSP < 50 GeV would 
already have been detected. Thus if the gaugino masses are given by anomaly mediation 
(with /i <C ttia), then we cannot rely on Wino pair annihilation, and instead we need a 
mechanism for suppressing C if this scenario is to be viable. However, because this result 
seems tantalizingly close to giving the correct relic density, we will consider below whether 
deviations from the basic anomaly mediated predictions are possible. 

3.3 Possibilities for Suppressing Moduli-Induced Gravitinos 
3.3.1 Approximate Symmetries 

If the modulus is charged under an approximate symmetry, then its decay rate to gravitino 
pairs may be naturally suppressed. For this to make sense, the modulus cannot have a shift 
symmetry, so it cannot be a volume modulus as in KKLT, but it may be a shape modulus. It is 
also important that there exists a sector of light particles with the same symmetry properties, 
so that modulus decay to this sector is unsuppressed - otherwise the overall modulus decay 
rate would be small, but the branching fraction to gravitinos would be unchanged. 
For instance, as a simple toy model we could consider 



= j d 4 e xtx + $t$ + $t$ c + e ^! + H.c.-^^ + . 

+ J d 2 B [A + fi 2 X + m^$ c ] (32) 



where $ and $ c have opposite charge under an approximate U(l) symmetry whose violation 
is parameterized by e. In the limit that e — > 0, the modulus cannot decay to gravitinos by 
charge conservation, so for small e the decay rate is suppressed. We could easily include a 
sector of light 'charged' particles so that the overall modulus decay rate is not suppressed. 
Whether or not such a situation could arise depends on the details of very high-energy physics, 
but there is no reason to expect that such a setup is impossible. 

3.3.2 Light Moduli 

If the modulus is light, with mass close to 2m3/2, then the branching fraction to gravitinos 
will be phase-space suppressed. Thus parametrically we would have 



lm\ - 4m 2 3/2 

C oc -5!—- -, (33) 

m<t> 

but obtaining C ~ 10~ 3 in this way would require a large, poorly motivated fine-tuning. 

A more interesting possibility is that there are two or more moduli with different masses. 
If the lightest modulus is lighter than 2m 3 / 2 , then it will not produce gravitinos, but its 
decay will dilute the gravitinos from earlier processes. If the two moduli have masses < 
2777.3/2 <C 777$, then the lighter modulus begins rolling later in the universe's history, so that 
p<j> ~ P* while they are both oscillating. After <3> decays, the lighter field dominates the 
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energy density of the universe until it in turn decays, releasing more entropy. Between the 
time of $ decay and the time of decay, radiation red-shifts as R~ A but the energy in only 
red-shifts as R~ 3 , so decay dilutes heavy relics such as gravitinos by a factor 

\ 3 / 4 /TV. \ 3 / 4 / \ 9 / 8 

(34) 



C ~ — ^ . (35) 



which could optimistically be as large as 1000. For instance, with m$ = 5 x 10 4 TeV, 
= 100 TeV, m 3 / 2 = 60 TeV, and N ~ 10, one would need C ~ 1/20 to obtain the correct 
relic abundance of dark matter (with these values direct production of dark matter from 
decay is not a problem). We will see that the presence of such a light modulus is also useful 
for alleviating cosmological bounds on the axion. 

Finally, we could simply return to the original scenario of Moroi and Randall [13] , with a 
single modulus with m<f> < 2m 3 / 2 . As they showed, the correct relic density can be obtained 
rather easily by relying on the pair annihilation studied in section 3.1 (in addition, they 
claimed that the modulus branching fraction to Winos was chirally suppressed, but this is 
false [33], [35]). It is certainly possible that moduli with masses of order m s simply do not 
exist, but this is a UV sensitive issue. 

3.3.3 Model Dependence 

Although models with C ~ 1 are generic, there are models where C is parametrically different. 
For instance, as shown in [35], in specific supergravity realizations of the KKLT scenario, one 
finds 

\ 77 

In this case ~ ~^~R as usua l ; but '0' is not a mass eigenstate. The mass eigenstate 
modulus mixes with the supersymmetry breaking sector, so that it effectively has F eigen ~ 

m 2 

—^-F giving a suppressed — > 2-0 3 / 2 decay rate. This result depends on specific properties 

of the supersymmetry breaking sector (the Polonyi field must be light, with mass of order 
7773/2) that need not be true in more general KKLT-type scenarios, but these theories are 
technically natural, and in some cases the terms that would give C ~ 1 can be forbidden by 
global symmetries. 

We have been assuming that the gauginos get there mass from anomaly mediation alone, 
but as mentioned above, it is possible that the same operators 

d 2 9-^-W a W a L D -^AA (36) 

M pl M pl v ; 

that allow moduli to decay to gauginos also contribute to the mass of the gauginos. If 
F^/Mpi ~ m^/rn^ as expected generically, then the Wino mass can be altered, and for a 
light gravitino and heavy modulus the Wino could become lighter. In that case, Wino pair 
annihilation could be sufficient to achieve the correct relic density; for instance, with 777^ = 70 
GeV and 7773/2 = 100 TeV, one obtains the correct abundance. Unfortunately this is a UV 
sensitive question, and it seems to require fine tuning to make the Wino light. Furthermore, 
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there is the possibility of introducing CP violation through these types of interactions. Note 
that these operators are negligable when is small and the <fi — > 2^3/2 decay rate is small, 
which is the domain of the section 3.3.1. 

A better solution is to use a large \x term in equation ( 1101) to make the Wino light. We 
would need to abandon the approximate PQ symmetry explaining the top-bottom hierarchy, 
but with n ~ vtla one obtains a lighter Wino, making it possible to explain the dark matter 
relic abundance purely in terms of pair annihilation. 

A final possibility is that supersymmetry breaking is dominated by D terms. As shown 
in the appendix, in this case C would be suppressed as 

Coc^V. (37) 



This occurs because only the sgoldstino decays to gravitino pairs, and in the limit of D- 
breaking, there is no sgoldstino (instead the goldstino has a vector partner, which is not 
of interest here). As shown in [37] it is possible to obtain parametrically large D terms, 
and Z)-breaking naturally gives a split supersymmetry type spectrum. This possibility may 
warrant further investigation, as it may help to alleviate the cosmological moduli/gravitino 
problem in more general circumstances. 

3.4 The 0^0 + ^3/2 Channel 

For completeness, we consider the decay channel <fi — > <f) + ip3/ 2 [11] . although we will see that 
it is dangerous, so we will need to forbid it kinematically. It comes from operators of the 
form 

L d [ d 4 6—X&$ 

M pl 



a- 



M 



P i 



- V3am 3 / 2 mi<jxj) + ■■■ (3£ 



where x is the goldstino, X = x + Q\ + 6 2 m^/2M p i is the field that dominantly breaks 
supersymmetry, and a is an 0(1) coupling constant that must be relatively large so that the 
modulino will be lighter than the modulus. As a complex scalar field, <fi actually has two real 
modes, with masses 

I \/3 

«V± = ym 2 ~± V3am 3/2 m^ w ± —am 3/2 (39) 

Generically, both modes will be present, but 0_ cannot decay to a modulino and a gravitino. 
If we require that is greater than the sum of the modulino and gravitino masses, then 
we must have a > (2 + m 3 / 2 / m ^)/ if this relation is violated, then the decay will be 
kinematically prohibited. 

With these assumptions, one can calculate the modulus decay rate and branching fraction 
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into modulinos and gravitinos [TT], obtaining 



v 2 ™^ ( _ K + m 3/2 ) 2 \ 3/2 / _ K-m 3 / 2 ) 2 \ 1/2 
^ m\ ) \ m l J 

^ (40) 

for the decay rate into a gravitino and a modulino, where in the second line we take 7773/2 <C 
77i0. Thus we find a relic abundance [11] 

We know that Ylsp ~ V 3 / 2 , so we obtain 

»w iS P > 10- cev x (sag-) (^i^y ( 1000 TeV ) 

LAP L6F ~ V100 GeV/VlOO TeV/ V "V / 



3/2 

(42) 



for the dark matter energy density of the universe. We expect that va^spYhSP = 4 x 10~ 10 
GeV [36] if the LSP accounts for all of the dark matter in the universe, so this decay mode will 
overproduce LSPs. Fortunately, we can prohibit this decay kinematically without introducing 
any fine-tuning. 



3.5 Weakening Cosmological Bounds on the Axion 

The axion [26] decay constant F is bounded from below due to astrophysical constraints, and 
in generic cosmological scenarios it is bounded from above by the requirement that the axion 
does not overdose the universe. The moduli fields in our model decay when the universe 
is at a temperature near Aqcd, releasing a significant amount of entropy, and potentially 
relaxing constraints on F by diluting the axions. Throughout we will be considering the relic 
axion energy density from misalignment production, which is by far the dominant method of 
production for large F. 

As shown in [27], diluting relic axions is a bit delicate. This is because there are three 
processes that need to be considered: a particle or field is decaying, releasing entropy and 
making the universe cool more slowly; H is decreasing, alleviating axion hubble friction; and 
the axion mass is increasing, since m a is strongly dependent on temperature because it arises 
from instanton effects. The axion only begins to roll as H drops below m a , so if our modulus 
decays before this point then the entropy released does not decrease the final axion density. 
In general, particles decaying after the universe has cooled below 1 MeV are dangerous to 
BBN, leaving a narrow window of 1 MeV < Tr < for reheating temperatures of decays 
that can dilute axions. Note that Tr should be interpreted as the approximate temperature 
of the universe after modulus decay - the temperature of the universe never increases. 

First we estimate Th, the reheating temperature above which the modulus does not dilute 
the axion density at all. As the modulus decay completes, the universe will be radiation 
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dominated with temperature Tr, so H ~ T^/M p i at this time. The axion has mass 

m a (0) = 13 MeV (^^) (43) 



at temperatures less than about Aqcd/k and 



m a (T) = 0.1m o (0) (^^^J . (44) 

at temperatures above Kqcd/^- If Tr = Th, then 3H(Th) = m a {Th) just as the modulus 
decays, so we find that 

T ft = 1.5 GeV f J (45) 

where we have taken Aqcd = 200 MeV. We see that modulus decay will certainly dilute 
axions with F near the current bound, but for F = 10 15 GeV, Th ~ Tr for the 1000 TeV 
modulus field. However, we can do much better with a lighter modulus. 

In order to obtain a relic axion energy density less than the current matter density in the 
universe, it was shown in 1271 that we must have 



f < 2xl0 »GeV(l^XV / 7^V /2 (46) 



T, 



R 



where T a is the temperature of the universe when m a (T) = 3H. We can estimate T a as 

r^OOMevf^rV^r (47) 
V 100MeV / V F ) 

thus as anticipated above, = 1000 TeV leads to no improvement. However, if we consider 
moduli with ~ 100 TeV as considered in section 3.1 and 3.3.2 then we obtain Tr ~ 6 
MeV, giving F < 7 x 10 14 GeV. 



4 Displaced Gluino Vertices at the LHC 

We know that cosmology favors a light anomaly mediated spectrum, so we can be justifiably 
optimistic that our model will be tested at the LHC. In fact, we expect that gluinos will 
be copiously produced, and that gluino pairs will decay at secondary vertices separated by 
distances of order a millimeter or more. To begin to analyze this process, note that the 
resolution of the LHC's ATLAS detector [S] will be 

73 

Ad = 1H »m (48) 

115 

A^o = 87 + 3 am (49) 

(pr/GeV)(sin(0))i 

Since gluinos are produced by QCD reactions such as gg — > gg, we can expect a large pt 
and an order one 9. Thus we will need gluino decay vertices displaced by > 10 fim to have a 
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chance of distinguishing them from primary vertices. There are four factors relevant to the 
displacement distance: the gluino production rate, the gluino lifetime, the relativistic time 
dilation, and the possibility of detecting gluinos from the tail of their distribution, which 
decays exponentially with distance from the primary vertex. 

We will not attempt to compute the the gluino production rate precisely [2T], [28], [29] . 
[30] , but we know that the answer is large, and our results will be very insensitive to the 
details. It was found in (9] that for m g ~ 350 GeV, the LHC will produce about one gluino 
per second, or about 3 x 10 7 per year of operation. For m g ~ 2 TeV, the LHC will produce 
at least a thousand per year. Thus we make the conservative assumption that the LHC will 
produce 

w =i*(i™y (so) 

V m 9 J 

gluinos per year. 

Various groups [5], [20] have examined the decay of the gluino and computed its lifetime 
in Split Supersymmetry, including the relevant one-loop operator running from the SUSY 
breaking scale to the TeV scale. The result is 

4 sec / m, \* (\ TeV x 5 



T « = — iwiw' I — I (51) 



m- g 



where K < 1 is a weakly varying function of m s , rrig, and tan(/3). This is extended by a 
factor of approximately 

3 TeV 

7 « (52) 



m-g 



due to relativistic time dilation. We translate the gluino lifetime into a displacement by 
multiplying by a factor of ca/1 — m 2 / (3 TeV) 2 ~ c. 

Putting these results together, we find that the number of gluinos with displacement in 
the interval (R, R + dR) is 



where the canonical distance Do is given by 



1 TeV \ 6 / m s 



Co=(7 . 2mm) ^_J [ _^_ ) (54 ) 

If we require the observation of at least a few displaced gluino vertices in one year's worth 
of LHC data, we find the limit 

i« 2 (r, „, /2TeV\\ /l TeV\ 6 / m s y 

10 < ( 7 + 4 ' 0g {—) ) {—) ( lOOOTw) (55) 

on the mass parameters of the model. This result is only accurate to within a factor of a 
few, but this is unimportant because displaced gluino vertices will be visible for virtually 
the entirety of our anomaly-mediated parameter space, as seen in figure 1. For comparison, 
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Gluino Mass (TeV) 

Figure 1: We will observe displaced gluino vertices at the LHC if m s and m g lie in the region above 
the curve. 



if we require the observation of 100 displaced gluino vertices each year at the LHC, then 
the digit 7 in the above equation changes into a 2, but this is the only alteration. Note 
that backgrounds for this signal from muons and b-quarks will be small due to the large jet 
energies, large missing energy, and especially because of the large vertex displacements. It 
was shown in the detailed study of [23] that simply using cuts on and missing Et, for 
vn-g < 1.4 TeV, the signal from R-Hadronized gluinos could be seen in 30 fb~ l of LHC data. 

5 Conclusions 

We have shown that a simple Split Supersymmetry spectrum based on anomaly mediation 
can satisfy cosmological constraints with a light spectrum. This is possible because the 
successive decays moduli — > LSPs and moduli — > gravitinos — > LSPs generate dark matter 
very efficiently, in contrast with most anthropically motived models [3], [18], [22] based 
on perturbative freeze out, which often require 1-2 TeV LSP masses ([23] is an exception). 
Unfortunately, in one of the most elegant scenarios dark matter is overproduced, but we found 
several mechanisms that can lead to the correct relic density, and the original mechanism of 
Randall and Moroi [13] remains a very viable possibility. We have also seen that our model 
will have a striking LHC signature - displaced gluino vertices. 

As outlined in the introduction, there are only a few clues for would-be anthropic model 
builders. The two most tantalizing are probably gauge coupling unification and dark matter, 
and we must assume that these two issues are resolved in concert if we are to avoid an almost 
infinite set of possibilities for new physics. Here we have shown that retracting the usual 
assumptions about the generation of dark matter do not necessarily make models irrelevant 
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for experimental collider physics, as we might have fearedZJ. Furthermore, we noticed a 
fortuitous accident - as a side effect, our mechanism weakens cosmological constraints on 
axions, expanding the rather narrow window on the axion decay constant. With so little 
information about new physics, we should be appreciative when a small piece of the strong 
CP problem falls into our lap. 

The success of anthropic predictions of the cosmological constant, combined with the 
derth of electric dipole moment signatures and new flavor changing neutral current effects, 
the smallness of the S and T parameters, and the lack of new physics at LEP point toward 
an anthropic solution to the two naturalness problems of high-energy physics and cosmology. 
If the world is supersymmetric at high energy, then the spectrum of new particles that we 
have considered is an excitingly predictive model for LHC physics. 



A Appendix: Moduli- Gravitino Decay Rates 

These results were obtained in collaboration with Aaron Pierce and Jesse Thaler. Similar 
methods can be found in [3Tj . 

Consider an arbitrary, globally supersymmetric lagrangian L((f>, if), F, A, A^, D, d^) that is 
a function of off-shell chiral multiplets and vector multiplets. Now assume that supersymme- 
try is spontaneously broken by some combination of D component and F component VEVs. 
The sgoldstindfl is the linear combination 

4>x = Yl "P"^' ( 56 ) 

i 

where F 2 = Eji^ 2 + S a D^ is the goldstino decay constant. Note that there is no reason to 
expect that <p x is a mass eigenstate. We will show that the only two derivative operator 
involving a scalar and two goldstinos is 

l^d^xa^d.x - jd v <j> x X^^- (57) 

Thus the rate of moduli decays to gravitino pairs is determined by the overlap of the mass- 
eigenstate moduli with <p x . As a simple corollary, we see that if supersymmetry were broken 
entirely by D terms, then there would not be a sgoldstino, so the rate for modulus decay to 
gravitino pairs would be zero. 

To prove this assertion, we need to do a field redefinition to introduce the non-linearly 
realized goldstino x- As in the conceptually simpler case of a goldstone boson, we can do this 
by choosing a 'vacuum alignment', and then parameterizing a supersymmetry transformation 
with the goldstino x( x )- For instance, if the theory consisted of a single chiral multiplet 
$ = (0, ip, F), we would represent $ = e l<5x (0, 0, F). 

7 However, it is worth noting that when dark matter is not generated by perturbative freeze out, we lose 
the elegant parametric prediction that the weak scale is the geometric mean of the Planck and Cosmological 
Constant scales [2J. The numerical prediction for the value of the cosmological constant is unchanged. 

8 If there are non-vanishing D component VEVs then there will also be a vector partner of the goldstino, 
but since we are only interested in moduli-goldstino interactions, we are ignoring these fields. 
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In the general case, when supersymmetry is linearly realized the goldstino is the combi- 
nation 

^ = E^ + E^- ( 58 ) 

i a 

To isolate it, we begin by rotating the fermions ipi, \ a into a mass eigenstate basis, so that 
the fermion lagrangian becomes 

L f = i^d^ g + Y, [i^d^ a + ^(M a + VM)] + - (59) 

a 

where the elipsis denotes interaction terms. Note that the superpartners of the fermions ip a 
have no auxiliary component VEVs. 

Now we perform a x parameterized supersymmetry transformation on L, and then choose 
the 'vacuum alignment' condition ip g = 0. We are only interested in two operators, the 
interaction term mentioned above and the kinetic term for x, which we require for canonical 
normalization. Both of these operators involve two goldstinos, so they can only come from 
C5 X A or C5 X A5 X B where A, B, and C are some combination of fields in the lagrangian. But 
the first type of term would be X3<^A and this vanishes on-shell because it is proportional 
to the A equation of motion. Thus only terms of the second kind are relevant. 

The supersymmetry transformation rules are the familiar 

5 X <P = V2 X *P (60) 

S x ij = iV2^x ] D^ + V2xF (61) 

5 X F = iy/2x^D^ (62) 

5 X \ = a^xF^ + ixD (63) 

5 X A^ = -i\^a^x + iX^A (64) 

6 X D = -x^D^ -D^ X ] - (65) 

where is the appropriate gauge covariant derivative. The x kinetic term can only come 
from the transformation of the ip g kinetic term, since the transformation of all other opera- 
tors involve too many fields. Furthermore, it is remarkable that two derivative interactions 
between a scalar and two goldstinos also only come from the ip g kinetic term, because the 
other terms either give too many fields or vanish on the equations of motion. Thus as claimed 
we find 

L D i X ] ^d, X + ^d^x^d.x ~ jd v <t> x X ] ^° v X ] + - (66) 

after canonically normalizing Xi where the elipsis denotes other interactions, including a 
two goldstino interaction with the vector superpartner of the goldstino. Although from the 
transformation rules it might seem that x has non-derivative interactions, these all cancel - 
X has a shift symmetry, just like a goldstone boson. Thus yukawa couplings with x cannot 
be generated. 

We should note that this is an effective field theory, and it is only valid for E 2 F = 
m 3 / 2 M p i. In most examples of interest (KKLT moduli, or moduli that get masses from 
supersymmetry breaking), the modulus mass easily satisfies this criterion. In any case, above 
this scale supersymmetry breaking is a small effect, and in place of a single goldstino one 
would need to consider the detailed dynamics of the hidden sector. 
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